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MAGNETIC FIELD ANALYSIS METHOD AND PROGRAMS 
FOR ROTATING MACHINES 



BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to a method of 
analyzing the magnetic field of a rotating machine. 

5 Description of the Related Art 

As described, for example, in a paper 
"Induction motor analysis by time-stepping techniques", 
T. W. Preston, A. B. J. Reece and P. S. Sangha, IEEE 
Trans, on Magnetics, vol. 24, No. 1, pp. 471-474, 1988, 

10 a conventional rotating machine magnetic field analysis 
method adopts a time-stepping method of sequentially 
analyzing the magnetic field of a rotating machine by 
stepwise rotating the rotor. 

In analyzing the magnetic field of a rotating 

15 machine, the rotor is stepwise rotated and a matrix 
equation having potentials as unknown variables is 
solved by a numerical solution approach such as a 
finite element method, generally by using an iterative 
solution method.- Since a rotating machine is 

20 accompanied by a magnetic saturation phenomenon, a 

permeability changes as the function of a magnetic flux 
density so that iterative calculations are required 
which are inherent to non-linear analysis. With the 
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above-cited rotating machine magnetic field analysis 
method, each time the rotor is stepwise rotated, a 
solution is calculated by the iterative solution method 
by using the solution obtained at the preceding time- 
step as an initial value of an unknown variable. In 
obtaining a solution with this method, however, it is 
necessary to suppress a rotation angle width narrow to 
some extent. This method becomes ineffective when the 
rotation angle width exceeds a certain value. To avoid 
this, the calculation is generally made by setting the 
initial value to 0. This requires, however, a large 
number of iterations necessary for obtaining a 
solution, taking a large calculation time. 

SUMMARY OF THE INVENTION 

An object of this invention is to provide a 
rotating machine magnetic field analysis method and its 
program capable of shortening a calculation time. 

According to one feature of a rotating 
machine magnetic field analysis method of this 
invention, a magnetic field distribution in the whole 
analysis space is obtained by using, as initial values, 
a value of a magnetic field distribution in the stator 
space and a value of the magnetic field distribution in 
the rotor space. 

Other features of the invention are given as 
stated in the appended claims. 

Other objects, features and advantages of the 
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invention will become apparent from the following 
description of the embodiments of the invention taken 
in conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
5 Fig. 1 is a diagram showing an example of a 

flow illustrating a rotating machine magnetic field 
analysis method according to the first embodiment of 
the invention. 

Fig. 2 is a diagram showing an example of a 
10 potential distribution on a slide plane. 

Fig. 3 is a diagram showing a potential 
distribution on a slide plane of components separated 
into each mode. 

Fig. 4 is a diagram, showing an example of a 
15 flow illustrating a rotating machine magnetic field 

analysis method according to the second embodiment of 
the invention. 

Fig. 5 is a diagram showing an example of a 
computer system. 

20 Fig. 6 is a diagram showing an example of a 

magnetic disk. 

Fig. 7 is a diagram showing an example of a 

CD-ROM. 

Fig. 8 is a diagram showing a typical example 
25 of displaying a convergence state of magnetic field 
'analysis by the invention. 



DESCRIPTION OF THE EMBODIMENTS 

Embodiments of a magnetic field analyzing 
method of this invention will be described with 
reference to the accompanying drawings. 

Fig. 1 illustrates an example of a flow of a 
magnetic field analysis method according to a first 
embodiment of the invention. Since a magnetic field 
distribution is analyzed generally by utilizing a 
potential, the embodiments will be described also by 
utilizing a potential. 

At the first step 11, a potential on a slide 
plane between the rotor and stator is calculated from a 
potential distribution at a predetermined time, e.g., 
n-th time t = nAt where At is a time-step width. This 
solution is separated into each mode along a rotation 
direction to derive constant components and rotation 
fundamental mode components. 

Fig. 2 shows an example of a potential 
distribution on the slide plane, and Fig. 3 is a graph 
showing the potential distribution separated into each 
mode. As shown in Fig. 3, the potential solution on 
the slide plane can be resolved into constant 
components, rotation fundamental mode components and 
harmonics components. As shown in Fig. 3, the constant 
components and rotation fundamental mode components are 
a main part of the potential components describing the 
rotation magnetic field. 

At the second step 12, an analysis space is 
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separated into a rotor space inclusive of the rotor and 
a stator space inclusive of the stator. 

The third step is divided into two steps, a 
3. 1-th step 13 and a 3. 2-th step 14, either one of 
which may be executed first. At the 3. 1-th step 13 a 
fundamental mode on the slide plane is rotated by a 
rotation angle of a rotation magnetic field 
corresponding to a time-step width and the constant 
components are added to the rotated fundamental mode 
components. This addition result is used as the 
boundary conditions to perform a non-linear magnetic 
field analysis by taking into consideration the 
magnetic saturation in the stator space. In this case, 
the solution obtained at the preceding time t = nAt is 
used as the initial value of a permeability 
distribution. 

At the 3. 2-th step 14, the rotation 
fundamental mode on the slide plane is rotated by an 
angle obtained by subtracting a rotation angle of the 
rotor from the rotation angle of the rotation magnetic 
field corresponding to the time-step width and the 
constant components are added to the rotated 
fundamental mode components. This addition result is 
used as the boundary conditions to perform a magnetic 
field analysis by taking into consideration the 
magnetic saturation in the rotor space. Also in this 
case, the solution obtained at the preceding time t = 
nAt is used as the initial value of a permeability 



distribution. In a synchronous machine, the rotation 
speed of a rotation magnetic field is equal to the 
rotation speed of the rotor. Therefore, the boundary 
field on the slide plane changes only finely because of 
the influence of harmonics components so that the 3.2- 
th step 14 may be omitted. 

At the 3. 1-th step 13 and 3. 2-th step 14, the 
solution at the preceding time t = nAt may be used as 
the initial value of an unknown variable or another 
initial value such as 0 may also be used- One example 
of the unknown variable is a magnetic vector potential 
A defined as B = rot A. In the finite element method, 
a projected integer value aj of A upon a mesh side is 
an unknown numerical value in terms of numerical 
analysis. The potential A can be expressed by A = 
EajNj (Nj is a vector base function). 

At the fourth step 15, by using as the 
initial values the potential solutions obtained at the 
3. 1-th step 13 and 3. 2-th step 14, a magnetic field 
analysis is again performed in the whole space by 
taking the magnetic saturation into consideration. 

In the non-linear magnetic analysis taking 
the magnetic saturation into consideration, the nearer 
the initial value of the potential is to the solution, 
the smaller the number of iterations for the iterative 
solution is and the shorter the calculation time is. 
In a synchronous machine among others, the 3.2.-th step 
14 can be omitted so that the calculation of obtaining 



the initial value near to the solution is only the 
analysis at the 3. 1-th step in the stator space so that 
a suitable initial value can be obtained with a smaller 
scale of calculations. 

Fig. 4 illustrates an example of a flow of a 
magnetic field analysis method according to a second 
embodiment of the invention. The second embodiment 
will also be described by utilizing a potential. 

At the first step 11, a potential on the 
slide plane between the rotor and stator is calculated 
from a potential distribution at an n-th time t = nAt 
where At is a time-step width. From this potential on 
the slide plane between the rotor and the stator, 
constant components and rotation fundamental mode 
components are derived. 

At the second step 12, the analysis space is 
separated into the rotor space inclusive of the rotor 
and the stator space inclusive of the stator. 

The third step is divided into two steps, a 
3. 1-th step 13 and a 3. 2-th step 14, either one of 
which may be executed first. At the 3. 1-th step 13 a 
rotation fundamental mode on the slide plane is rotated 
by a small angle such as 1/2 or 1/3 the rotation angle 
of the rotation magnetic field corresponding to a time- 
step width. The rotated fundamental mode components 
added with the constant components are used as the 
boundary conditions on the slide plane to perform a 
non-linear magnetic field analysis by taking into 



consideration the magnetic saturation in the stator 
space. In this case, the solution obtained at the 
preceding time t = nAt is used as the initial value of 
a permeability distribution. 

At the 3. 2-th step 14, the rotation 
fundamental mode on the slide plane is rotated by a 
small angle such as 1/2 or 1/3 the angle obtained by 
subtracting a rotation angle of the rotor from the 
rotation angle of the rotation magnetic field 
corresponding to the time-step width. The rotated 
fundamental mode components added with the constant 
components are used as the boundary conditions on the" 
slide plane to perform a magnetic field analysis by 
taking into consideration the magnetic saturation in 
the rotor space. Also in this case, the solution 
obtained at the preceding time t = nAt is used as the 
initial value of a permeability distribution. The 3.2- 
th step 14 may be omitted from the reason stated 
earlier . 

At the 3. 1-th step 13 and 3. 2-th step 14, the 
solution at the preceding time t = nAt may be used as 
the initial value of an unknown variable or another 
initial value such as 0 may also be used. In order to 
obtain a solution in a short time, it is preferable to 
use the solution at the preceding time t = nAt. 

At the fourth step 15, from a change in the 
potential solutions- and permeability distributions 
obtained at the 3. 1-th step 13 and 3. 2-th step 14, the 



potential analysis and permeability distribution at the 
time t = (n+l)At can be estimated through linear 
extrapolation. By using this solution as the initial 
value, a magnetic analysis is again performed in the 
whole space by taking the magnetic saturation into 
consideration . 

In the second embodiment, since the rotation 
angle of the fundamental mode is smaller than that of 
the first embodiment, a change amount of the solution 
and permeability distribution at the time t = nAt is 
small. Therefore, the calculation time at the 3. 1-th 
step 13 and 3. 2-th step 14 is relatively short. It is 
possible to obtain an approximate solution and 
approximate permeability distribution at the time t = 
(n+l)At at a higher speed than that of the first 
embodiment . 

In the embodiment, an approximate value at 
the time t = (n-fl)At is obtained through linear 
extrapolation based upon one analysis in the rotor and 
stator spaces. Instead, the approximate value at the 
time t = (n+l)At may be obtained by analyzing the two 
cases in the rotor and stator spaces through two- 
dimensional function extrapolation when the fundamental 
mode is rotated by two different rotation angles 
smaller than the rotation angle width corresponding to 
the time-step width At. In this case, although the 
calculation amount is doubled, the approximate value 
becomes nearer to the correct solution so that the 
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whole analysis at the time t = (n+l)At can be obtained 
in a shorter calculation time. 

In the magnetic field analysis, a convergence 
solution of a large scale coarse matrix equation is 
5 generally obtained by an iterative solution method such 
as an ICCG method (a conjugate gradient method with 
incomplete Cholesky decomposition) . Also in this 
embodiment, the solution is obtained by the iterative 
solution method. In both the embodiments, as shown in 

10 Fig. 8, the convergence state of calculations can be 

checked by calculating the equation with a computer 80 
and displaying the calculation result on a display 81 
connected to the computer 80. The display contents in 
Fig. 8 show a gradual change in remaining differences 

15 of the magnetic field analysis solutions in the stator 
space shown in Figs. 1 and 4 obtained by the iterative 
solution method. The remaining difference is defined, 
for example, by |Ax - b | / |b | for a matrix equation of 
Ax = b. 

20 Although only the potential constant 

components and rotation fundamental mode components are 
used, the harmonics components may also be used in the 
similar manner. In this case, the approximate solution 
can be made more precise so that a solution by the 

25 whole space analysis can be obtained in a shorter 
calculation time. 

According to the magnetic field analysis 
method for a rotation machine described above, the 



analysis is made by dividing the analysis space into 
the rotor and stator spaces so that an approximate 
solution can be obtained at high speed. Since the 
whole space analysis is made by using solutions nearer 
to the correct solution, the final solution can be 
effectively obtained faster than a conventional method. 
This effects become larger if the analysis system is 
large . 

The embodiments described above may be 
realized by a dedicated system. As "illustratively 
shown in Fig. 5, the embodiments may be realized by a 
general computer system and a program running on this 
system, the computer system having: a keyboard 51; a 
computer 52 constituted of input means for inputting 
data described earlier and the program, a storage unit 
for storing the input data and program, a calculation 
unit and the like; and a display 53. 

When the program is to be supplied to such a 
computer system, the program is stored in a storage 
medium such as a magnetic disk 61 such as shown in Fig. 
6 and a CD-ROM 71 such as shown in Fig. 7. The storage 
medium distributed and kept is installed in the 
computer 52 to read the program with a magnetic disk 
drive or a CD-ROM drive of the computer 52. If the 
program distributed via a communication network is 
input to the input means, the program is stored in the 
storage medium such as a magnetic disk to allow 
repetitive use of the program. 
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According to the invention, it is possible to 
provide a magnetic field analysis method for a rotating 
machine and its program capable of shortening a 
calculation time. 

It should be further understood by those 
skilled in the art that although the foregoing 
description has been made on embodiments of the 
invention, the invention is not limited thereto and 
various changes and modifications may be made without 
departing from the spirit of the invention and the 
scope of the appended claims. 



